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Abstract 

An algorithm is given to compute a normal form for hyperelliptic 
curves. The elliptic case has been treated in 0]. In this paper the hyper- 
elliptic case is treated. 

1 Introduction 

The problem we consider is the same as in || , but then for the hyperelliptic case 
instead of the elliptic case. Let C be an algebraic curve given by an irreducible 
polynomial / G C[x,y]. The function field of C is the field of fractions of 
C[x,y]/(f), and is denoted as C(x,y). We want to decide if C is hyperelliptic 
or not. If C is hyperelliptic, we want to find a Weierstrass normal form for C 
(i.e. an equation yg 2 = F(xq) for some squarefree polynomial F(xq) G C[xo]). 
We want to find a birational map from C to this normal form, and its inverse. 
This means expressing x, y in terms of Xq, yo and vice versa. Let g be the genus 
of C. 

The problems to be solved are the following 

• Decide if C is hyperelliptic or not. By definition, C is hyperelliptic if g > 1 
and there exists Xq G C(x,y) such that C(xq) is a subfield of C(a;,y) with 
index 2. If C is hyperelliptic then the field K := C(xo) is unique. 

• Find a generator xq of K. 

• Find an element yo of C(x,y), not in K. Take for example yo — x (inter- 
change x, y if it turns out that x € K). Compute the minimal polynomial 
of yo — x over K. This is done by writing xq as a quotient V\/v2 with 
v\, V2 G C[x, y] and computing the resultant ReSj^i^io — vi> f) where xq 
is the name of a variable that plays the role of Xq. Then add an element 
of K to yo to change its minimal polynomial into the form y^ — F(xq) for 
some rational function F(x ) G K. Then multiply y by an element of K 
to turn Fq into a squarefree polynomial. 



• At this point the functions xq, yo are expressed as elements of C(x, y) and 
their algebraic relation j/g — F(xq) is in Wcierstrass form. Now view xq, yo 
as variables xo, yo, and express x, y as rational functions in these variables. 
This can be done in the same way as in M . 



2 Decide if C is hyperelliptic 

If / S C[x, y] is irreducible, then we can decide if C is elliptic by computing the 
genus g. The curve is elliptic if and only if g = 1. For the hyperelliptic case 
more work needs to be done. An algorithm for this is implemented in Maple 6. 
To view the it, issue the following commands in Maple 6 or higher 
interface (verboseproc=2) : op( 1 algcurves/is_hyperelliptic ' ) ; 
and the following code appears 

proc(f, x, y) 

local a, b, i, j , g, z; 

option remember, 'Copyright (c) 1999 Waterloo Maple Inc. All rights\ 
reserved. Author: M. van Hoeij 1 ; 
i := degree (f, {x, y}) ; 

j := min(degree(f , x) , degree(f, y) , i - ldegree(f, {x, y})); 
if i < 4 or j < 2 then false 

elif j = 2 then evalb(l < ' algcurves/genus ' (f , x, y)) 
else 

b := 'algcurves/diff erentials' (f , x, y, skip_dx) ; 
g := nops(b) ; 

j := l/2*(i - l)*(i - 2) - g; 

if g = 2 then true 

elif g < 2 or j < 2 or j < 4 and 

{seq(z[3], z = ' algcurves/singularities 1 (f , x, y))} = {1} 

then false 

else 

z := {seq(a[i] , i = 1 .. g)}; 

j := {coeff s(expand(numer(subs(RootOf (f , y) = y, evala( 
Expand ( 

subs(y = RootOf(f, y) , add(a[i] *b[i] , i = 1 .. g)~2)))))) 

, z)>; 

z := {seq(a[i] , i = 1 .. nops(j))}; 

z := nops(j) - nops(indets (subs ('solve/linear 1 ({coeff s( 

expand (add (a [i] *j [i] , i = 1 .. nops(j))), {x, y})>, z) 
, z)) intersect z) ; 
if z < 2*g - 1 then error "should not happen" 
else evalb(z = 2*g - 1) 
end if 
end if 
end if 
end proc 
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We will not explain all steps, but the main steps are the following. 

• Let V be the vector space of holomorphic differentials. An algorithm 
described in Q is used to compute a basis b±, . . . , b g for V. 

• If g < 2 then the curve is not hyperelliptic and the algorithm stops. 

• Let W be the span of all bibj, 1 < i < j < g. 

• If C is hyperelliptic, and if t is any generator of K, i.e. K = C(t), then it 
is known that 

fdt, i = 0,... , 5 -l (1) 

is a basis of V. Then i*d 2 t, % = 0, . . . ,2g — 2 must be a basis of W. So if 
C is hyperelliptic then the dimension of W must be 2g — l. It is not hard 
to show that if C is not hyperelliptic and g > 1 then dim(VT) > 2g — 1. 
In the algorithm, z = dim(W). Testing if z = 2g — 1 decides if C is 
hyperelliptic or not. 

3 Compute a function of degree 2 

Once we have a generator xq for K, it is not difficult to find yo, see section [I]. 
This is implemented in Maple in the procedure 1 algcurves/genus2 ' (to view 
the code type similar commands as in the previous section). After that, the 
remaining steps are the same as in Q . 

So the main remaining problem to be solved is to compute some generator xo 
of K. We have a basis b\, . . . ,b g oiV. We want to find two particular elements 
of V, whose quotient is a generator of K. For example, if t is any generator of 
K, then dt, tdt will do, and so will t 9 ~ 2 dt, t 9 ~ 1 dt. We will describe how to find 
a subspace V2 C V of the form Ct 9 ~ 2 dt + Ci ff_1 dt for some generator t of K. If 
i>i, t>2 is any basis of V%, then we can take xq = V1/V2 as our generator of K. 

The idea for finding V% is the following. Take an arbitrary nonsingular point 
P on the curve. Since t may be any generator of K , we may assume without 
loss of generality that t(P) = 0. Then t 5_1 df has the highest valuation at P of 
any nonzero element of V, and t 9 ~ 2 dt has the second highest valuation. So we 
can find V% from V by solving linear equations: As long as dim(y) > 2, let m be 
the smallest valuation at P of any element of V. Then compute the subspace of 
V of those elements that have valuation > m, and replace V by that subspace. 
This way the dimension decreases, and after a number of steps V will equal V%. 
Then take xq — V\ /V2 where V\ , V2 is any basis of V% . This is implemented in 
Maple in the procedure ' algcurves/index2subf ield' . 

4 An example 

f = y 9 + 3x 2 y 6 + 3x 4 y 3 + x e + y 2 
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This curve has genus 3 and is hyperelliptic. The result of the Maple commands: 

with(algcurves) : Weierstrassf orm(f ,x,y,xO,yO) ; 

is 



[y 2 + 2 - 7x + 35 x 4 - 21 x Q 5 + 7 x 6 - x 7 + 21 x 2 - 35 x 3 , 

(-1 + x ) 2/o, 1 - 3 x + 3 x 2 - x 3 



y(y 3 + x 2 + l) x (y 3 + x 2 ) 3 



y 6 + 2 x 2 y 3 + x A + y ' jy 

The first expression is the Weierstrassform y 2 — F(xo). The second expression 
is a generator for the field K, i.e. a function of degree 2 in C(x,y). Under the 
birational map this corresponds to xq- The third expression is the image of yo 
under this birational map. The 4'th and 5'th expression are x, y expressed in 
terms of xq, yo- 
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